Spiral Waves in Media with Complex Excitable Dynamics 



Andrei Goryachev and Raymond Kapral 
Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, ON M5S 3H6, Canada 

The structure of spiral waves is investigated in super-excitable reaction-diffusion systems where 
the local dynamics exhibits multi-looped phase space trajectories. It is shown that such systems 
support stable spiral waves with broken rotational symmetry and complex temporal dynamics. The 
main structural features of such waves, synchronization defect lines, are demonstrated to be similar 
to those of spiral waves in systems with complex-oscillatory dynamics. 



Studies of spatially-distributed active media have 
demonstrated the ubiquity of self-organized spatio- 
temporal patterns, in particular spiral waves, in vari- 
ous physical or biological systems such as the Belousov- 
Zhabotinsky (BZ) reaction 0], catalytic surfaces 0), car- 
diac muscle |§ and colonies of the amoebae Dictiostelium 
discoideum [01. 

Most research has been devoted to the study of simple 
oscillatory or excitable systems. Recently it was shown 
that reactive media with complex periodic and chaotic os- 
cillations are capable of supporting spiral waves || with 
a variety of distinctive features, absent in simple oscil- 
latory systems ||. The rotational symmetry of spiral 
waves in period-doubled media is broken by synchroniza- 
tion line defects where the phase of the local oscillation 
changes by multiples of 2tt. It was conjectured that spiral 
waves with broken rotational symmetry could also be ob- 
served in super-excitable systems where the phase space 
trajectory after excitation follows a multi-looped path of 
relaxation to the stable fixed point Q. Broken spirals 
with a clearly visible synchronization defect line ema- 
nating from the spiral core were observed under special 
three-dimensional conditions in the BZ reactive medium 
. The nature of spiral waves in complex-excitable me- 
dia is nevertheless largely unexplored. 

In this paper we show that spiral waves with broken ro- 
tational symmetry exist in a prototypical super-excitable 
system and demonstrate that the topological properties 
of line defects, described earlier for oscillatory media, also 
hold for excitable systems. 

We consider a spatially-distributed system whose dy- 
namics is governed by a pair of reaction-diffusion equa- 
tions of the form 



du 1 

-wr = -"I 1 

at e 

dv 

m =u - v 



u)(u 



-f(v)) + DV 2 u, (1) 



where v is a non-diffusive variable and both u and v are 
functions of time and space. This model with f(v) = 
was studied in (9) as a simplified version of the FitzHugh- 
Nagumo model which serves as a prototype of an ex- 
citable system described by two variables. The excitable 
dynamics of system (Q) consists of two fast and two 
slow stages. If displaced from the stable fixed point 
u = 0, v — to the right of the unstable branch of the 
nullclinc u = 0, it quickly reaches upper stable branch 



u = 1. It follows this branch until v reaches sufficiently 
large values and then jumps to the lower stable branch 
u = 0, along which it slowly relaxes to the stationary 
state. To add a super-excitability to ([!]) a modification 
of the unstable branch of the nullcline ii — was pro- 
posed in (if]] as 
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The introduction of (g) changes the shape of the unstable 
branch of u = so that for suitably chosen parameters 
a, a and Vq it nearly touches the nullcline v — (see 
Fig.|l|). As a result, if another excitation is applied to 
the system (Q) before it has reached the stationary state, 
it may execute second, smaller excitable loop before it 
finally reaches the stable state. 
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FIG. 1. Phase plane (u,v) plot of a 1 type trajectory cal- 
culated at a point in the medium shown in Fig.^. Two stable 
branches (u = 0, u = 1) and one unstable branch of the null- 
cline u = are shown as solid lines and the nullcline v = by 
a dashed line. 

The spatio-temporal dynamics of a one-dimensional ar- 
ray of such elements forced by an external pacemaker 
with varying period was studied in |frO| . When the pe- 
riod of forcing Tf is larger than a certain internal period 
of the system To, the response is a train of waves corre- 
sponding to the large relaxation loop (1°). If Tf < Tq 
the system develops wavetrains with low amplitude cor- 
responding to small relaxation loop (0 1 ). Non-trivial be- 
havior is observed when Tf is only slightly smaller than 
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To . In this case the system shows mixed-mode waveforms 
(l n ) consisting of one large and n small waves. Response 
of this type is an example of the complex-excitable dy- 
namics targeted in our studies. Instead of an external 
pacemaker we use a spiral wave as a self-sustained source 
of excitation in the medium. 




FIG. 2. An irregular spiral wave with complex excitable 
dynamics is shown at two different times in the upper and 
lower panels, respectively. The concentration field v(r,to) is 
color-coded as shown in FigQ 

Spiral waves were initiated in a two-dimensional square 
domain with no- flux boundary conditions. While a and 
a were used as bifurcation parameters, other parameters 
were fixed at (e = 0.005, a = 0.6, b = 0.03, v = 0.2). 
As in jl0| complex-excitable dynamics was found in the 
parameter region between domains of large-amplitude 1° 
(small a and a) and small-amplitude 1 (large a and a) 
waves. In this region the medium supports mostly ape- 
riodic, stable spiral waves lacking rotational symmetry. 
Figure || shows such a spiral wave for a = 0.15, a 2 = 



0.001 at two time instances. Note that the wave length 
of a large-amplitude wave is larger than that of a small- 
amplitude wave and, thus, the shape of the spiral is dis- 
torted. The concentration time series v(r,t) at different 
locations in the medium (cf FigJ^) show aperiodic con- 
catenations of l 1 and l 2 waveforms, while trivial pat- 
terns 1° and 1 are completely absent. 




FIG. 3. Time series v(r, t) calculated at a point in the 
medium shown in Fig.^. 

As the domain of complex-excitable dynamics is tra- 
versed from 1° to 1 the contribution of the waveforms 
l n with n > steadily grows, as well as the number n 
of the small-amplitude loops. Thus, the transition from 
large to small-amplitude spiral waves occurs gradually 
through a succession of irregular spiral patterns with 
a progressively growing contribution of low-amplitude 
waves. 

Although irregular patterns are found in most of the 
complex-excitable domain, pure period-3 dynamics was 
found in a sub-domain of this region. Figure ^ shows 
the spatial structure of a spiral wave in this parameter 
region. Observation of the spiral wave dynamics for long 
time periods shows that the entire concentration field 
slowly rotates around the spiral core with a constant an- 
gular velocity u). The period-3 dynamics is manifested in 
a coordinate frame centered at the spiral core and rotat- 
ing with velocity to. Indeed, it takes three rotations of 
the spiral for the concentration field to return to itself. 

Figure |I shows a phase portrait of the dynamics at a 
non-special location in the medium calculated in the ro- 
tating frame. Before closing onto itself the phase space 
trajectory executes two small loops and one large loop, 
corresponding to a pure l 2 dynamics. Consider a po- 
lar coordinate frame (p, 4>) in the (u, v) phase plane with 
origin at an arbitrary point internal to both the small 
and large loops. During one full period of the dynam- 
ics the phase variable (f> changes by 6tt. Calculation of 
the phase at every point in the medium at a time to gives 
an instantaneous snapshot <fi(r,to) of the time-dependent 
phase field <f>{r,t). 
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FIG. 4. Spiral wave with period-3 dynamics. Color-coding 
is the same as in Fig.^. 

Consider a closed contour F that surrounds the spiral 
core. The phase increment Ar</> = / r V<f>(r, to)dl along 
r will be equal to a multiple of the full 6ir period of the 
dynamics. From the topological theory of point defects 
jnj in simple oscillatory or excitable media it follows 
that Ar</> is invariant and for the one-armed spiral waves 
in this study takes the values ±2n. In complex-periodic 
media this contradiction is resolved by the existence of 
synchronization defect lines emanating from the core || . 
The phase of the local oscillation experiences jumps equal 
to multiples of 2ir when such line defects are crossed. 
Any contour T encircling the spiral core intersects these 
lines so that the total phase increment Ar</> is obtained 
from the integration of V</>(r,to) along T yielding ±2tt 
plus phase jumps at the intersections of T with the syn- 
chronization defect lines. The sum of both contributions 
yields the full period phase increment of the local dynam- 
ics. It was predicted || that one should be able to find 
phenomena analogous to synchronization defect lines in 
media with complex-excitable dynamics. 

For the specific case of period-3 complex-excitable dy- 
namics where the total phase increment is 6n, one might 
expect to find one line defect where the phase jumps by 
An or two defect lines where the phase jumps by 2tt on 
crossing each line. In Fig.|| one sees two lines emanating 
from the spiral core at an angle of ~ 180°. Investiga- 
tion of the change in phase of the local dynamics across 
these lines shows that they are indeed 2ir synchronization 
defect lines exhibiting the loop exchange phenomenon de- 
scribed earlier for complex-oscillatory media. 

Figure ^ shows the v(r, t) time series at four neighbor- 
ing locations in the medium along a path traversing one 
of the defect lines. Panel a is a plot of the normal l 2 
dynamics seen on one side of the defect line. Every third 
minimum is lower than the two preceding minima and 
corresponds to the larger relaxation loop in the phase 
space plot. 
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FIG. 5. Change in character of local dynamics across one 
of the defect lines seen in Fig.^. See text for details. 

As one approaches the defect line, the large-amplitude 
loop shrinks while both small-amplitude loops grow 
(Fig.||(b)). Then, one small-amplitude loop begins to 
grow faster than the other and at some point be- 
comes larger than the still shrinking large-amplitude loop 
(Fig.||(c)). This loop exchange process continues until 
the new large loop attains a size equal to that of a large- 
amplitude loop and the other two loops shrink to the 
size of the small-amplitude loop. The local dynamics on 
opposite sides of the defect line (compare Figs ||(a) and 
|](c)) experiences a 2tt phase shift. Thus, the total phase 
increment Ar4> — 2n resulting from the integration of 
V<£(r, to) along T, excluding its intersections with defect 
lines, plus the two 2n phase jumps at the intersection 
points gives the expected 6tt phase increment. 

Our results show that the structure of complex- 
periodic spiral waves is governed by general topological 
principles independent of whether the dynamics is ex- 
citable or oscillatory. This fact allows one to extend 
the predictions inferred from the studies of complex- 
oscillatory systems to systems with super-excitable dy- 
namics. The formation of complex-periodic spiral waves 
in such systems might play a role in the development 
of some pathological conditions in the heart. Indeed, 
mixed-mode electrical activity with alternating large and 
small amplitude maxima, so-called alternans, is typically 
observed as a symptom of tachycardia [Q . 
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